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Abstract 

The Lund area law was developed into a Monte Carlo 
program LUARLW. The important ingredients of this gen- 
erator was described. It was found that the LUARLW sim- 
ulations are in good agreement with the BEPC/BES R scan 
data between 2-5 GeV. 

1 INTRODUCTION 

The hadron production mechanism in particles collisions 
is one of the important subjects in the study of strong inter- 
action. Quantum chromodynamics (QCD) is considered as 
the theory of strong interaction. However, the hadroniza- 
tion processes belong to nonperturbative problem for which 
no practicable calculation based on the first principle avail- 
able. Some phenomenological hadronization models were 
thus built up, which play important roles in the studies to- 
ward the final understanding of strong interaction. The fa- 
mous Lund string fragmentation model is one of the suc- 
cessful hadronization schemes, which contains several non- 
trivial dynamical features and describes the general semi- 
classical picture of hadron production. At high energies, 
the Lund generator, JETSET, can simulate the processes of 
hadron production via single photon annihilation and pre- 
dicts the many properties of the final states correctly. But 
the application of Lund model at intermediate energies has 
been blank. A direct way out of this situation is to start 
from the basic assumptions of Lund model and find the so- 
lutions of the area law without adopting any high-energy 
approximation. Based on the Lund area law, a new gener- 
ator LUARLW was compiled, which agrees with BES data 
between 2 — 5 GeV well (see Figure |l] on page 

2 LUND STRING FRAGMENTATION 

The foundations of Lund model (relativity, causality and 
quantum mechanics) are universal. The basic hadron pro- 
duction picture is string fragmentation. The produced new 
pairs (qq) and (qqqq) may form mesons and baryons if they 
carry with the correct flavor quantum numbers, otherwise 
they just behave like the vacuum fluctuations and do not 
lead any observable effects in experiments (see Figure 
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Figure 2: String fragmentation by a set of new pairs (qq) 
and [qqqq] production, hadrons form at vertices 

Using the assumptions of very high energy approximation 
(the remaining string always has large energy scale), left- 
right symmetry (fragmentation from go end or go end are 
identical) and iterative fragmentation (string fragmentation 
may be treated iteratively), Lund fragmentation function 
f{z) was derived uniquely, 

/(.)^^(i-.roxp(-^), (1) 

where, a and h are fragmentation parameters, m is the 
(transverse) mass of fragmentation hadron, and z is the 
fraction of light-cone momentum. J{z) is used in JET- 
SET to govern string fragmentation. Lund fragmentation 
function /(z) has the characteristics of inclusive distribu- 
tion, and the single particle production is independent of 
anything else before and after. The applicable region of 
f{z) is the remnant string still has large invariant mass. At 
intermediate energies, the mass-shell conditions should be 
the component part of the fragmentation dynamics, and the 
string usually fragments into 2^6 hadrons. Therefore 
the string fragmentation have to be treated as exclusive one 
instead of inclusive like in JETSET. 

3 LUND AREA LAW 

Lund string fragmentation process is Lorentz invariant 
and factorizable. The finite energy (s) system containing n 




Figure 3: The situation after n steps fragmentation. 

hadrons may be viewed as a cluster of infinite string frag- 
mentation system with energy (sq —> oo) (see Figure ||). 
According to the general properties of iterative cascade, 
the combined distribution is the product of fragmentation 
functions for every steps 

n 

dpn = Y[fjizj)dzj 

J=l 

= Cn ■ dpcxt{s,z) ■ dpn{ui,- ■ ■ ,Un), (2) 

where, C„ is normaUzation constant. We know from ^ 
that a subsystem may be split up from the total system, the 
processes occurring in the subsystem is the same as it be 
a complete system starting at the some original energy s. 
The external part 

dz 

rfpcxt (s, z) = rfs— (1 - z)" • exp(-6r) (3) 

z 

corresponds to the probability that the cluster will occur. 
The internal part 

/ 

\ ^ = 1 

n 

• Yl d'^PojS{p% - mlj) cxp(-6^„) (4) 

then corresponds to the exclusive probability that the clus- 
ter will decay into the particular channel containing the 
given n particles with energy-momentum {poj} and noth- 
ing else. -4icst — -4„ is the area enclosed by the quark 
and antiquark light-cone energy-momentum lines of n par- 
ticles. The factor 



|M|2 = exp(-5A0 



(5) 



may be viewed as the squared matrix element, the other 
parts are phase-space elements. In formulas, b is fun- 
damental color-dynamical parameters. Distribution M) is 



Figure 4: The vertex V divides the n-body string fragmen- 
tation into two clusters which contain ni and n2 hadrons 
and with squared invariant masses si and §2 separately. 



called Lund area law. The total area 

-4tot 



J2A,^An + r, T=-!^L_J1 (6) 



and 




(7) 



Finishing the integral over kinematic variables of 
n— particles, area law has following forms: 



• String ^ 2 hadrons 

C2 



P2 



• String ^ 3 hadrons 



exp(-6^^^^) + exp(-M^^^) 
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String ^ 4, 5, 6 hadrons 



dsids2 



: exp(-6r) 



\/X{s,Sl,S2) 
■ pni{si,Al)pn2{s2,A2)- (10) 

In above fragmentation distributions, the gluon effects 
are neglected. At intermediate and low energies, the emit- 
ted gluons from initial quark or antiquark are usually soft, 
most of which will stop before the string starts to break, 
the effect of the gluon will then essentially be small trans- 
verse broadening of two-jet system, the gluon and quark 
will then look as single quark jet. The gluon emissions do 
not significantly change the topological shapes (sphericity 
and thrust) of final states, and therefore no observable jet 
effects. 



4 MULTIPLICITY 



6.1 Scheme I 



Lund area law may give the expression of the muhiplic- 
ity distribution of fragmental hadrons. Define dimension- 
less n-particle partition function 

Zn = s J dRn ■ exp(-6^), (11) 

where dRn is the n-particle phase space element. The re- 
lation between Z„ and the multiplicity distribution P„ for 
primary hadrons is 



p — 

Pn has the approximative expression 

Pn = — j- • exp[co + ci(n - ^) 



(12) 



+C2{n-ix)\ (C2<0). (13) 



Quantity /i may written as the energy-dependent form phe- 
nomenologicaly 



or 



fi = a + h- exp(c\/s), 



H = a + 61n(s) -|- cln^(s). 



(14) 



(15) 



All parameters a, b, c, co, ci, and C2 need to be determined 
by experimental data and have been tunned with BES data 
samples of R scan. 

5 EXCLUSIVE DISTRIBUTION 

There are some different production channels for n- 
particle states, such as 4-body states may be 7r"'"7r~7r"'"7r~, 
7r"'"7r~7r°7r°, p'^p~TT + 7r~, etc. The exclusive probability 
for the special channel is 



In Lund model, quantum mechanical tunneling effect is 
was used to explain the production of new pairs qiQi. Parti- 
cles obtain their transverse momenta from the constituents. 
At each production point the (gjgi)-pair is given ±q^ and 
the particle momenta are 

P_L1 = qi, • • -P-Lj =<lj - qj-l, • ■■P±n = -fin- (17) 

Based on the Lund model, the following distribution with 
forward-backward symmetric correlation was derived 



F(")(qi,...,q„) 



: Cn exp 
: Cn exp { - 



1 

1 

4^ 



'2 , (qa - P2qi)^ 

qi "T 



I -Pi 



2 I 2 
qi + q„ 



with 
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(18) 



The correlations pj are phenomenological parameters, 
which in general are small. The covariant matrices 
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Vn 
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P2---Pn-1 P3---Pn-1 
\ P2---Pn P3--- Pn 



P2P3 ■■■ Pn \ 
P3--- Pn 



Pn = Bn- (VPS) ■ (SUD) ■ p„(m_Li, . . . , m±n; s). (16) give the correlations (pxiPxj) and (pyiPyj). 



• Bn is the combinatorial number stemming from may 
be more than one string configurations lead to this 
state. 

• (VPS) is the vector to pseudoscalar rate. 

• (SUD) is the strange to up and down quark pair prob- 
abiUty. 

6 TRANSVERSE MOMENTUM 
DISTRIBUTION 

Above results are obtained when the transverse momen- 
tums of all primary hadrons have given. In LUARLW, two 
transverse momentum distributions were used alternatively. 



6.2 Scheme II 

An available Gaussian-like transverse momentums dis- 
tribution as options in LUARLW for n— particles fragmen- 
tation reads 




(19) 



The conditional distribution of p_Li, p_L2, . . . , p_Lj are 




A^i • exp 
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(n - i + 1) (P-Lj + E; 



_P±i_\2' 
n=l n-j+l' 



{n - 3) 



2^2 



The final p_L„ is determined by energy-momentum conser- 
vation.The effective variance of p_Lj is 



The transverse momentum distribution in this scheme is not 
exact Gaussian type due to the transverse momentum con- 
servation and the threshold conditions. 



The well-know Monte Carlo simulation packet JETSET 
is not built in order to describe few-body states at the few 
GeV level in e+e" annihilation as at BEPC. We develop 
the formalism to use the basic Lund Model area law di- 
rectly for Monte Carlo program LUARLW, which will be 
satisfied to treat two-body up to six-body states. In LU- 
ARLW, the effects of all gluonic emissions were neglected. 
The LUARLW predicts more than 14 distributions totally 
agree with BES data well. 
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7 SUMMARY 
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Figure 1: e+e hadrons spectrum of raw BES data (hatched region) and LUARLW/SOBDRUNK (black line) at 
Ecm = 2.2 GeV. 



